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Abstract. Localized scattering phenomena may result in the formation of 
stationary matter waves originating from a compact region in physical space. 
Mathematically, such waves are advantageously expressed in terms of quantum 
sources that are introduced into the Schrodinger equation. The source formalism 
yields direct access to the scattering wave function, particle distribution, and 
total current. As an example, we study emission from three-dimensional Gaussian 
sources into a homogeneous force field. This model describes the behaviour of an 
atom laser supplied by an ideal Bose-Einstein condensate under the influence of 
gravity. We predict a strong dependence of the beam profile on the condensate size 
and the presence of interference phenomena recently observed in photodetachment 
experiments. 
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1. Introduction 

Sources of particles are known to be an indispensable element in all scattering 
experiments. Located far away from the scattering region, sources will give rise to 
boundary conditions on the scattering wave function in the form of incoming plane 
waves. However, sources that are close to or within the relevant interaction region 
should be introduced directly into the Schrodinger equation [1]. More generally, the 
concept of a particle source arises naturally if a complex quantum process can be 
decomposed into several steps, where the preparation of the state under consideration 
is separated from its further evolution in an external field. For example, in intense-field 
laser-atom physics, the laser interacts with a pool of atoms in a complicated nonlinear 
way. As a result of that process a coherent source of above-threshold electrons is 
generated that will emit electrons in the form of an outgoing spherical wave with a 
given energy in the continuum [2, 3, 4]. 

In this article we discuss in more detail a stationary quantum source that allows 
for a fully analytical solution. Our novel example predicts how the properties of the 
resulting matter waves are controlled by the size of the underlying source. The model 
we are presenting is surprisingly simple: A weak, harmonic perturbation couples the 
bound state in an isotropic oscillator potential to a continuum state experiencing a 
homogeneous force field. In section 2, we show how the ground state wave function 
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acts as a stationary quantum source for the uniformly accelerated scattering wave, 
and discuss the ensuing integrated current and current distribution. For a three- 
dimensional Gaussian source embedded into a linear potential environment, an exact 
solution for the emitted wave and the associated current is obtained in section 3, 
which allows for a pictorial interpretation in terms of a "virtual" point source. In 
section 4, this theory serves as a simple model for two phenomena recently investigated 
experimentally, the formation of a freely falling "atom laser" beam from a trapped 
Bosc-Einstcin condensate, as well as near-threshold photodetachmcnt microscopy of 
negative ions in an electric field. Finally, a brief outlook is given in section 5. 



2. Quantum sources in an external field 

To motivate the quantum source approach, we will introduce a two state model with 
time-independent Hamiltonians. One state is a bomid state of the Hamiltonian Htrapj 
the second state is a scattering state of the Hamiltonian i?cont containing an external 
potential. Both states have an energy difference of AE = Scont — -2'trap and are weakly 
coupled by a homogeneous but oscillating interaction potential of strength hQ: 

{indt - i?co„t)^co„t(r, t) = ftfle-'^^'/V^raplr, t) (1) 

{ihdc - Ftrap)VA,rap(r, t) = ftfie+'^^*/ Vcont(r, t) . (2) 

We split off the time dependence of the states 

V'co„t(r, t) = e-^^--*/Vcont(r) (3) 
V-traplr, t) = e-^^--*/ Vtrap(r), (4) 

to obtain the stationary equations 

(-Econt - -f^cont)V'cont(T') = Sl^V'trap ) (5) 
(£^trap - fftrap)V'trap(r) = hnipcont{r). (6) 

Upon introducing a suitable energy-dependent Green function Gcont('"5 r'; E) for i?cont 
(see section 2.2) 

{E - ffcont)Gcont(r, r'; E) = 5{r - r'), (7) 
the formal solution of equation (5) is given by 



t(r) = h^j dV C?cont(r, r'; E,,^,) V'trap(r'). (8) 

A similar equation holds for tptra-pir)- However, if we assume only a weak interaction 

in the sense that ),''t,rap('") is not changed appreciably by the interaction, we may as a 
first approximation replace V'trap('") in equation (8) by the bound eigenstate ipoir) of 
Htrap, which is defined by 

(£o - -fftrap)V'o(r) = 0. (9) 
In this way we have decoupled both equations and we obtain a new Schrodinger 
equation with an inhomogeneous source term a{r) = hfltpoir) 

(-Econt - i?cont)'0cont(T') = (^{r)- (10) 

Prom now on we will suppress the label cont, since we are not interested in the 
propagation of the bound state. Therefore the rewritten equation (8) 

'tP{r) = J d^r'G{r,r';E)a{r'). (11) 

serves as the starting point for our development of the theory of quantum sources. 
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2.1. Currents from sources 

We associate a probability current density j(r) with the scattering wave i'{r) via 

ji^r) = -Immr)* VV(r)} - ^^\Mr)\\ (12) 
m m 

where A{r) denotes the vector potential. Together with equation (10) it is 

straightforward to derive the equation of continuity for this stationary problem: 

V-j{r) = ~lm{a{rriPir)}. (13) 

Note that the introduction of sources a{r) in the Schrodinger equation (10) causes the 

appearance of a source term in the equation of continuity that depends on the wave 
function ^(r). By integration over a surface enclosing the source we obtain a bilinear 
expression for the total probability current carried by 'tp{r): 

J{E) = -|lm jy dV y dV a(r)* G{r, r'; E) (T(r')| . (14) 

This quantity details the total cross section of the underlying scattering process. 



2.2. The energy-dependent Green function 

For our stationary problem, the mathematical formulation of source theory heavily 
relies on the energy-dependent Green function G{r,r';E) of the Hamiltonian. In 
the continuous spectrum of H, however, this integral kernel is by no means unique. 
From a physical point of view, this ambiguity of the solution set (11) is required 
to accommodate different boundary conditions for the resulting wave function 
Here, we are interested in outgoing matter waves that enforce the use of the retarded 
energy Green function [5] . (This choice of kernel guarantees the correct sign for the 
total current J{E) (14), as shown by an eigenfunction expansion of G(r, r'; E) [6].) For 
reasons of causality, the retarded energy Green function is connected to the propagator 
K{r,t\r',0) of the equivalent time-dependent problem by the Laplace transform: 

1 r°° 

G{r,r';E) = - dte'^'/'' K{r,t\r',0). (15) 
1" Jo 

The time-dependent quantum propagator is a thoroughly studied subject, and a 

fairly exhaustive list of available solutions is given in [7, 8]. In contrast, only few 
closed analytic solutions are known for energy-dependent Green functions in three 
dimensions. Notable examples are the expressions for the free field environment, the 
static uniform electric field [9, 10, 11], the static uniform magnetic field [12, 13, 14], 
and combined parallel static electric and magnetic fields [15, 16]. 



2.3. A sum rule for the total current 

Exploiting the time-reversal symmetry relation for the propagator K {r , t\r' , 0)''' = 
K{r,—t\r',0), the total current J{E) (14) can be rewritten using equation (15) as 

J{E) = ^ J°" dte'^*/^ J dV J <i^r'a{r)* K{r,t\r',0)a{r'). (16) 

Integration with respect to the energy E, together with the initial condition 
K{r, 0|r', 0) = 6{r — r'), leads to the following sum rule for the total current: 

J°° dEJ{E) = ^ Jd^rJ d^r'a{r)*K{r,0\r',0)a{r') = ^ J d^r\a{r)\^. (17) 
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3. Sources in a uniform force field 

To illustrate the formalism developed in the preceding section, we discuss the specific 
example of a quantum source of Gaussian shape embedded in a uniform force 
field F. Experimentally, this situation is approximately realized in near-threshold 
photodetachment of negative ions in an electric field [17, 18, 19, 20, 21], where the limit 
of point-like sources applies. Considerably extended sources emerge when ultra cold 
Rubidium atoms are continuously released from a trapped Bose-Einstein condensate 
under the influence of gravity ("atom lasers", [22, 23]). 

3.1. Propagator and the energy Green function 

The Hamiltonian ffcont for uniformly accelerated motion with A{r) = o is given by: 

Hfieid = ^ - r • F. (18) 
2m 

For simplicity, wc align the field along the 2:-axis: F = F Cz- The propagator assigned 
to this Hamiltonian is well known and reads [7, 8, 24]: 



Kfi,M{r,t\r',0) = (-—) exp<^- 



\27rifii/ h 



TTl . f .1) Ft , / > F t 

— \r- r'r H (z + z') 

2t, ' ^ 2 ^ ^ 24m 



(19) 



The energy-dependent Green function is also available (see [9, 10, 11] for a compact 
derivation) and may be expressed in terms of Airy functions [25] : 

777 1 

Gfieid(r, r'; E) = ^y^r^^ [Ci(a+)Ai'(a_) - Ci'(a+)Ai(a_)] , (20) 

where a± = -(][2E + F{z + z')±F\r-r'\], (3 = [m./ [Ah^F^)]'"'^ , and C\{x) = 
Bi(a;) -|- i Ai(a;). For further reference, we note that this Green function is invariant 
with respect to simultaneous shifts of the origin and the energy: 

Gfieid(r, r'; E) = Gfieid(r - r', o;E + Fz'). (21) 

For a compact notation it is suitable to introduce a set of scaled parameters 

^ = PFx p = jSFr 

V = (3Fy e = -2PE (22) 
C = /3Fz T = t/{2h(3). 
This allows us to express the Green function in integral form via (15) and (19): 

Gfieid(p, o; e) = -2\p{(3Ff — exp i^^p'' + iT(C - e) - ^ j • (23) 



3.2. Point-like sources 



Initially, wc discuss a point-like quantum source in a homogeneous force field F. 
Assuming an isotropic emission pattern, we may idealize the source cj{r) in terms of 
the Dirac ^-distribution 

cjs{r) = C5{r), (24) 

where C denotes some complex constant. We note that for this type of sources the 
sum rule stated in equation (17) is not applicable, since the £^ norm of 5{r) is not 
defined. Instead, the detailed mechanism behind the source in equation (24) merely 



Matter waves from quantum sources in a force field 5 

influences the emission process through the scaUng parameter C, the source strength. 
The wave function generated by (Ts{r) directly follows from equations (11) and (20) 

7P{r) = CGReidir,o;E) (25) 

and yields analytic expressions for the current density (12) in the field direction 

Ur,E) = \Cf£^ {z[A^ia.)f+(3[z{2E + Fz) + Fr'] [Ai{a.)f} , (26) 
and for the total current (14) 

J{E) = lim Im {Gfieid(r, o; E)} 

= ^^^J^ {[Ai'{-2(3E)f + 2(iE[Ai{-2(3E)f} . (27) 
In a different context, these expressions are implicitly contained also in [10, 15]. 

3.3. Gaussian source 

Now we turn our attention to more realistic spatially extended sources. Specifically, 
we consider a source term derived from the wave function V'o(^) in equation (9) of 
isotropic Gaussian shape 

£7(r) = mi/joir) = hriNo exp(-rV(2a2)). (28) 

The parameter a describes the width of the source and A^o = (i~^/^7r~^/^ denotes the 
proper normalization from the condition 



/ 



dV|Vo(r)r = l. (29) 



To obtain the expressions for the currents generated by a Gaussian source, we start 
with the derivation of the wave function from equation (11). Working in the time- 
dependent propagator representation (see equation (15)), we can carry out the r' 
integrations of Gaussian type in 

V'(r) = -iQNo r dte'^'/" [ d^r' KiiM{r,t\r',0)e-^"^^'^'''\ (30) 



Upon introducing the scaled width a = fiFa and shifted parameters 

C,=Q + 2a^ e = e + Aa^, (31) 
the problem is reduced to a single complex integration 

oxp -p^ + iuiC - e~) - ^ . (32) 

-2iQ2 (i7r?i) ' 12 / 

Here, A(g) = hn{2y/^af^^e^°'''^^-'^"'/^\ andp^ = ^2^,^2^^2_ We moved the temporal 
integration into the complex plane by substituting u = t — 2ao? . The representation 
chosen in equation (32) emphasizes the close relationship of ?/'('") to the Green fimction 
Gfieid(Pj o; e) in the form (23). To evaluate this integral analytically, we split the path 
of integration into two sections, one along the real u-axis, the other one along the 
imaginary u-axis: il>{r) = ipnearir) + ipfarir). The contribution due to 

Vnear(r) = -2 iA(g)/3(/3F)^ ^!L^e"^V«+i«(C-?)-i«Vi2 

„A(gW!v^,-,V,..", (33) 
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Laser wave number [cm Detuning Frequency ii = E/{27Th) [kHz] 



Figure 1. Characteristics of matter waves from sources in uniform force fields. 
Left panel: Total current J(E) in near-threshold photodctachment of S" in a 
homogeneous electric field F = 2.205 ■ 10*eV/m. (O ) Experimental data of 

Gibson et al [21]. { ) Theoretical result from (27) with adjusted source 

strength C. Right panel: Number of atoms remaining in a Bosc-Einstein 
condensate N{T) after continuous release of atoms for T = 20 ms as a function 
of the detuning frequency u. (O ) Measurement by Bloch et al reported in [26]. 

( ) Theoretical prediction according to (37) and (41) with effective Gaussian 

condensate width a = 2.8 /xm and outcoupling strength Q = 2it ■ 105.585 Hz. 



is a purely real term which drops off very quickly with increasing distance p from 
the source region. The more interesting far- field contribution ipfur^r) can be exactly 
evaluated using the integral representation (23): 

/•OO 1 

r^Ur) = -2iA{i)Pif3Ff / _^^^ipVu+MC-e)-iuVi2 

Jo (Ittu) ' 

= A(e)Gfieid(p,o;e). (34) 
With the help of equation (21) we can cast the last equation into the form 

^Ur) = hn{2V^af\-"^^'^/^'+^-'^'^y(^^'^ Gfieid(r, -^a'; E). (35) 

This expression displays a remarkable feature of the wave function ■(/'far('") originating 
from a Gaussian source: The extended Gaussian source can be replaced by a virtual 
point source of the same energy at a location shifted by mFa^/ (2?i^) from the centre 
of the Gaussian distribution, carrying the energy dependent weight A(e). From this 
relation, the expression for the currents due to (28) are conveniently found from 
the analogous expressions for a point source by just performing the indicated shifts. 
Neglecting ■i/'ncar('"), the far-field current density reads according to equation (26) 

X |J {c[Ai'(e~ - C + P)f + P [C(C - e) + p'] [Ai(e~ - C + p)f} . (36) 

The same procedure yields the total current. However, since both ipnca.T{f) and a{r) 
are purely real and only the imaginary part of ipir) is needed for the evaluation of 
the total current (see equation (13)), the following expression obtained by shifting the 
energy in equation (27) is even an exact result: 

J(e) = 64^'''/2fir!2a3/5e4"'(-~-4"V3) |[Ai'(?)]' - e~[Ai(e~)]2} . (37) 
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Figure 2. Interference fringes in uniformly accelerated motion. 
Distribution of electrons in near-threshold photodetachment of in a 

homogeneous electric field F = 423eV/m. Distance from the source z = 0.514m, 
electron energy E = 100.5 /ieV. Left panel: Image recorded by Blondel et al [20]. 
Centre panel: Calculation from (26). Right panel: Corresponding current density 
profile as a function of the lateral distance. 



In the limit of extended Gaussian sources, a simple approximation to this formula 
is available that leads to a geometrical interpretation. We start out with a time 
dependent integral formulation of the total current that follows from equations (14) 
and (23) after the spatial integrations are performed: 



J 



(in) 



3/2 



The integral is evaluated in saddle point approximation. Assuming a ^ e, we keep 
only the leading order terms of a Taylor expansion in the exponent and prefactor. The 
resulting current has Gaussian shape: 

2^^M^^^_,./(4„.)^ (39) 
a 

This expression is equivalent to the implicit representation 

Jsp{E) = ^J<i^r \a{r)\^ 6{E + Fz). (40) 

Evidently, the approximation (40) obeys the sum rule (17) for the total current J{E). 
For extended sources, the energy dependence of Jsp{E) reflects the source structure: 
By the resonance condition E + Fz = 0, the total current probes the density |V'o(^')| 
on different slices across the source. 



4. Examples: Photodetachment and atom laser 

As an application for the developed formalism, we want to discuss two physical systems 
that have recently undergone experimental evaluation. 



4-1- Photodetachment 

In photodetachment experiments, negatively charged ions are illuminated by a laser 
beam to release the surplus electron. If the photon energy closely matches the electron 
affinity of the ionic species, the detached electron starts out with a well-defined 



Matter waves from quantum sources in a force field 



8 



minuscule amount of kinetic energy E. Then, the size of the emitting ion is small 

compared to the initial electronic wavelength A — h/\/2mE, and details of the atomic 
structure cannot be resolved: The detached electron wave is solely characterized by 
its orbital angular momentum (which is in turn fixed by the selection rules of the 
underlying dipole transition), and allows for a description in terms of a point source 
C7(r). For isotropic (s-wave) emission, the Dirac ^-distribution as{r) (24) is a viable 
choice. 

In a field- free environment, the electron is emitted simply into a spherical wave, 
tp{r) (X e''°''/r. In this case, the quantum source model presented in section 2 
recovers the well-established behaviour of the photodetachment cross section near 
threshold known as Wigner's law [27] which states that J{E) cx ^/E. Recently, the 
question how the threshold behaviour of the photocurrent is altered in the presence 
of external fields in the interaction region has found considerable attention [17, 28], 
and precision experiments were conducted to measure both the energy dependence of 
the cross section [18, 21] and the corresponding photoelectron distribution [19, 20] for 
photodetachment in a homogeneous electric field environment. Here, the electric force 
F = —eE accelerating the detached electrons allows for direct comparison with the 
results presented in section 3.2. The left panel in figure 1 displays the photocurrent 
measured by Gibson et al [21] together with the source-theoretical prediction (27) first 
obtained in a different fashion by Fabrikant [15]. The agreement is striking. Besides a 
slow onset below threshold indicating tunneling emission, the plot prominently features 
a "staircase" appearance of J{E) modifying Wigner's law. Although an explanation 
of this electric field induced effect in terms of the closed free-falling orbit has been 
offered [29] , it is more accurately described as the remaining imprint of a remarkable 
interference pattern in the current density distribution jz{r,E) (26). 

In a recent series of experiments [19, 20], Blondel et al recorded the spatial 
distribution of the photocurrent on a distant detector plane perpendicular to 
the applied electric field. One of the impressive images obtained using the 
"photodetachment microscope" is depicted in figure 2 (left panel), revealing an 
interference ring pattern of macroscopic size. In the source formalism, the electronic 
distribution should reflect the local current density jz{r,E) pertaining to the Green 
function Gficid('", o, £■) (20) in a imiform force field. Indeed, an image calculated 
from (26) closely resembles the experimental result (centre panel in figure 2). For 
comparison also the corresponding lateral current profile is shown (right panel). 

The conspicuous fringes find a simple semiclassical interpretation in terms of two- 
path interference: As first pointed out by Demkov [30], within the sector of classically 
allowed motion always two ballistic trajectories will connect the origin with a given 
location on the detector. Hence, the uniform electric field effectively acts as a two- 
slit interferometer [11], providing a sensitive device for the precise determination of 
electron affinities [20] . Similar results are expected for point sources in parallel electric 
and magnetic fields where four-path interference takes place [16]. 

4-2. Atom laser 

In section 2 we introduced quantum sources using a coupled two state model. This 
model may serve as a simplified description of the atom laser, a coherent beam of 
particles released from a trapped Bose-Einstein condensate (BEG) [22, 23]. The two 
states are realized by different hyperfine Zeeman levels of Rb. The state Vtrap('") is 
trapped by a magnetic field, whereas the continuum state V'cont('") is not infiuenced by 




Figure 3. Transition from a point-like source to an extended source. 

Depicted is the total current J{E) from a Gaussian source of freely falling Rb 

atoms for different source widths a versus the detuning frequency u. ( ) Exact 

currents from (37). ( ) Slicing approximation according to (40). The coupling 

strength is H = 27r • 100 Hz. Beam profiles are shown in figure 4. 



the trap but subject to gravitational attraction. These states are coupled through an 
oscillating magnetic field of adjustable frequency with a specific coupling strength hil. 
Near its minimum, the potential for the trapped state is approximately harmonic. For 
non-interacting particles, the ground state of the condensate in this potential is given 
by a oscillator wave function. In our discussion, we maintain the Gaussian shape of 
the condensate wave function ipoir), but increase its width a in order to account for 
the repulsive atomic interaction. 

The ensuing quantum source a(r) continuously releases Rb atoms (mass m) 
with an initial energy E = 2-Khu determined by the frequency of the perturbation. 
These particles are subsequently accelerated by the gravitational force F = mg 
{g 9.81 m/s^). From section 3.3 it is possible to predict the efhciency of this process 
(the total particle current) as well as the spatial distribution of atoms (their current 




Figure 4. Size dependence of tlie particle distribution. 

Profiles of the current density 0,z,E = 27ri/?l) for Gaussian sources emitting 
free-falling Rb atoms. All results are given for i/ = 2.5 kHz and Q = 2tt ■ 100 Hz. 
There is rotational symmetry about the z-axis (sec figure 2). The profiles are 
calculated from (36) at a distance z = 1 mm from the source. The source width 
is denoted by a. See also figure 3. 



density profile). Experimental data concerning the number of remaining condensate 
atoms N{T) after T = 20 ms of atom laser operation is available [26]. Obviously, this 
number is related to the current by 

N{T) = N{0)exp[-J{E)T], (41) 

where N{0) denotes the initial number of atoms in the BEC and J{E) the total current. 
Besides the gravitational force, the relevant parameters entering the theoretical 
prediction for J{E) in equation (37) are the coupling strength hfl and the Gaussian 
width a. In figure 1 (right panel), the calculated number of remaining atoms (41) 
is compared to the experimental measurement by Bloch et al reported in [26]. The 
coupling frequency fl used for the calculation is fixed by the sum rule (17) applied to 
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the experimental data, and the effective width a = 2.8 /xm (that is actually largely 
governed by atomic repulsion in the BEC) is obtained from a fit. Contrary to 
the case of photodetachment, the current characteristics faithfully reproduces the 
Gaussian shape of the source, as stated by the approximation obtained by slicing 
the condensate at height z ~ E/F (40). However, according to the exact expression 
for J{E) (37), source theory predicts a dramatic change in behaviour of the total 
current for smaller condensate sizes as illustrated in figure 3. The plots show both 
the analytic solution J{E) and the corresponding approximation Jsp{E) for various 
source widths a. According to the sum rule (17), the area under all curves is the 
same. For a > 1 /itm, both models yield almost identical curves. However, for smaller 
condensates the differences become noticeable and for a < 0.4 /im, the slicing model 
fails completely: J{E) becomes asymmetric with respect to v = 0, and oscillations in 
the current are prominent. These features are familiar from the photocurrent discussed 
in the previous section, where we related them to two-path interference in a uniform 
field (figure 2). 

To examine if similar interference fringes occur also in an ideal atom laser beam, 
wc compute the current density profiles for the set of Gaussian sources displayed in 
figure 3 from (36). The results from this calculation are shown in figure 4 where we 
choose as initial energy E = 2nhi' with u = 2.5 kHz. The distance from the source is 
fixed at z = 1 mm. In this figure, a distinct ring pattern in the current density prevails 
for a < 0.4 /xm. The number of fringes diminishes with increasing source width, until 
for a > 0.8 /itm the current profile attains an increasingly narrow Gaussian shape. To 
interpret this behaviour, wc first note that for an extended source region, the simple 
concept of two interfering paths originating from a single point in space is not readily 
applicable. Recalling the particular property of a Gaussian source to act as a virtual 
point source shifted in space (35), we may recover the concept of interfering paths. 
However, the efi'ective initial kinetic energy decreases with growing source size (31) 
and becomes negative for E < mF^a^/{2h?), leading to a "virtual" tunneling source 
that emits a beam of Gaussian profile. Its properties are discussed in detail in [11]. 

5. Conclusion and outlook 

Quantum sources that are inserted into the Schrodinger equation provide a convenient 
and practical tool to assess and solve scattering problems in external potentials. 
Specifically, wc studied point-like and Gaussian quantum sources embedded within 
a three-dimensional homogeneous force field, and analytical solutions for the resulting 
particle distribution and total current were obtained. For point sources, the model 
is in excellent agreement with data gained from near-threshold photodetachment 
experiments conducted in the presence of an electric field. Available measurements on 
continuous atom laser beams are in reasonable accordance with the theory presented 
for an extended source. Dependent on the condensate size, we predict the appearance 
of a ring pattern related to two-path interference in the atom laser beam released from 
a single ideal Bose-Einstein condensate. 

A modification of the presented theory to comprise sources with definite angular 
momentum is feasible and will be the subject of a forthcoming publication. Physical 
applications of this extension include photodetachment experiments involving p-wave 
emission and the effects of vortices in an ideal Bose-Einstein condensate on the 
outcoupled atom laser beam. 
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